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The rise in the energy spectrum of the positron ratio, observed by the PAMELA satellite above 
10 GeV, and other cosmic ray measurements, have been interpreted as a possible signature of Dark 
Matter annihilation in the Galaxy. However, the large number of free parameters, and the large 
astrophysical uncertainties, make it difficult to draw conclusive statements about the viability of 
this scenario. Here, we perform a multi-wavelength, multi-messenger analysis, that combines in 
a consistent way the constraints arising from different astrophysical observations. We show that 
if standard assumptions are made for the distribution of Dark Matter (we build models on the 
recent Via Lactea II and Aquarius simulations) and the propagation of cosmic rays, current Dark 
Matter models cannot explain the observed positron flux without exceeding the observed fluxes 
of antiprotons or gamma-ray and radio photons. To visualize the multi-messenger constraints, we 
introduce "star plots" , a graphical method that shows in the same plot theoretical predictions and 
observational constraints for different messengers and wavelengths. 

PACS numbers: 95.35. +d, 98.35.Gi, 98.35. Jk 



I. INTRODUCTION 

Dark Matter (DM) annihilation or decay can in prin- 
ciple produce significant fluxes of positrons, antiprotons, 
photons, neutrinos and other secondary particles. Re- 
cently, the positron channel has received a lot of atten- 
tion, since the PAMELA collaboration has released the 
data relative to the positron fraction [T that exhibit a 
spectacular rise, which is in agreement with earlier re- 
sults from AMS-01 [2] and HEAT [3l IH [5] , and compat- 
ible with the claimed 300—800 GeV excess in the elec- 
tron plus positron spectrum measured by ATIC-2 balloon 
flights [6 . Lying in the energy range > 10 GeV, such an 
excess has prompted a large number of papers putting 
forward explanations that include DM annihilations or 
decays in the galaxy, and nearby astrophysical objects 
like pulsars. 

Here, we consider the DM annihilation hypothesis 
and perform a multi-messenger analysis in order to con- 
strain the properties of viable DM candidates. More 
specifically, we study positrons, antiprotons, 7-rays and 
synchrotron emission due to the propagation of elec- 
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trons and positrons in the Galactic magnetic field. The 
multi-messenger approach has already provided useful 
constraints on DM scenarios. For instance, the non- 
observation of an excess of cosmic-ray antiprotons up 
to - 100 GeV by the PAMELA satellite indicates 
that hadronic annihilation channels should be strongly 
suppressed. This motivates a simple distinction between 
DM candidates: leptophilic, i.e. that annihilate mainly 
into lepton pairs, and hadrophilic, whose annihilation 
final states are gauge bosons or quark pairs and in- 
duce non-negligible fluxes of both positrons and antipro- 
tons. Obviously, the former override p bounds, while 
the latter need to be rather heavy in order to suppress 
p fluxes below ~ 100 GeV, but predict inevitably larger 
fluxes at higher energies, that should be soon probed. 
Other interesting messengers for DM searches are high- 
energy neutrinos. In fact, neutrino observatories such as 
Super-Kamiokande [8J and IceCube [9] are able to detect 
upward-going muons produced in the interaction of high- 
energy neutrinos within the Earth interior. Therefore, 
once the DM profile is fixed, neutrino observations of the 
Galactic Centre region effectively constrain the proper- 
ties of dark matter, especially for multi-TeV candidates 
- see e.g. [101 E]. 

Alternative detection strategies include searches for 
new physics in accelerators, and direct detection experi- 
ments, searching for keV nuclei recoils due to dark mat- 
ter scattering off nuclei (see e.g. Refs. [12]). A convinc- 
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ing identification of DM can probably be achieved only 
through a combination of different detection strategies. 

Now, if one assumes that DM candidates are thermal 
relics from the early universe, the present annihilation 
cross section needs to be aannV ^ 0(10"^^) cm^s~^ in 
order to produce the observed relic abundance flcdmh^ ~ 
0.1. Such a thermal relic with a TeV mass needs a 
~ 10'^ boost in the annihilation flux to accommodate 
the PAMELA excess T^, which can hardly be pro- 
vided by the "dumpiness" of the Galactic halo (see be- 
low for a thorough discussion). Non-thermal relics, non- 
standard cosmologies, or velocity dependent ("Sommer- 
feld enhanced" ) annihilation cross-sections, have been in- 
voked to circumvent this problem, but the question re- 
mains of whether the large cross-sections needed to ex- 
plain the positron data can be made consistent with other 
astrophysical observations. We perform here an exten- 
sive analysis of the multi-messenger constraints in the 
framework of the latest high-resolution numerical simu- 
lations of a Milky Way like halo. We explore two spe- 
cific classes of models: leptohilic candidates, inspired in 
Refs. [Ml [15], and hadrophilic candidates, with a specific 
emphasis on models of Ref. [TS]. For all these models 
we compute the flux of: i) positrons, ii) anti-protons. 
Hi) gamma-rays from the Galactic Centre, iv) gamma- 
rays from the Galactic halo, and v) synchrotron emission 
due to the propagation of electrons and positrons at the 
Galactic Centre. We take into account the dependence of 
the annihilation cross section on the relative velocity. In 
particular, we calculate the boost factor due to Sommer- 
feld enhanced substructures, in the framework of the Via 
Lactea II |17j and Aquarius [181 simulations, and we dis- 
cuss the consequences for the gamma-ray flux from DM 
annihilations in Galactic and extragalactic halos. 

The paper is organised as follows. In Sec. [IT] we spec- 
ify our prescriptions for the distribution of DM, includ- 
ing the mass-function and concentration of DM clumps. 
In Sec. |III| we introduce the DM candidates used in the 
analysis. Sections |I V| to [VTI] arc then devoted to compute 
fluxes in the different channels under scrutiny: positrons, 
antiprotons, 7-rays and synchrotron. The main conclu- 
sions are drawn in the last section. 



II. ASTROPHYSICAL INPUT 

Over the last few years, N-body simulations have im- 
proved considerably, and recently two groups have pub- 
lished the results of high-resolution simulations. Via 
Lactea II (VL2) [17 and Aquarius (Aq) [TH]. In the 
former, both smooth and clumpy components are well 
fitted by Navarro-Frenk- White (NFW) profiles and the 
abundance of subhalos follows the rather steep behaviour 
M~^, while in the latter the density profiles seem to be 
Einasto-like and a shallower subhalo abundance cx Af^^'^ 
is found. A common feature is the presence of many re- 
solved subhalos and a characteristic dependence of their 



concentration on the position inside the halo, as we will 
soon explain. 

The smooth density profile is well modelled, in the VL2 
and Aquarius scenarios, by: 
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being r the distance to the galactic centre. The local 
density is Pq = Psmifo) with = 8 kpc. Following 
[T^ . the NFW profile of Via Lactea II is considered to 
saturate at a radius Tgat such that Psmi'Tsat) = Psat = 
2 • lOi* Mokpc"^ 

The density profiles inside clumps are NFW in Via 
Lactea II and Einasto with a — 0.17 in Aquarius. The 
corresponding concentration parameters are fitted by a 
double power law in mass 
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where M is the mass of the clump and r again the galac- 
tocentric distance. For NFW and Einasto profiles, one 
has Ts - r2oo/c2oo and = M200/ (/o''°° dr' 4^r'V(^')) 
being r' the radial coordinate inside the clump, f{r') 
the radial dependence of the clump density profile and 
r2oo the radius which encloses an average density equal 
to 200 times the critical density of the universe. In this 
way, the clump inner density profile pci{M,r,r') is un- 
ambiguously defined once the clump mass M and the 
distance to the galactic centre r are specified. A quan- 
tity that turns out to be relevant for galactic positrons 
and antiprotons is the so-called annihilation volume, 
^(M,r) = /;-° dr'ATTr'\p,i{M,ry)/pQf. ^ 

Another important input from N-body simulations is 
the spatial and mass distribution of clumps: 



VL2 



Aq 



dMdV 



2 ' 



exp 



dMdV 

2 I I r 

Ra 



{M,r)^Ash{M/MQ)-'-^x 



- 1 



,a = 0.678 



in units of Mg^kpc^^. In the following, we will refer to 
these expressions as Psh- 

The normalization Ash is fixed according to the find- 
ings of numerical simulations. In the VL2 case, we im- 
pose that 10 % of the galaxy mass M^r is virialized in 
structures with mass in the range [10~^M„jr, 10~^M„ir]. 
In the case of Aquarius, we require that 13.2 % of Myir is 
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FIG. 1: (Color online) The quantity Nci{£,)ai as a function of 
Mmin for Via Lactea II (solid red line) and Aquarius (dashed 
blue line). The vertical line indicates the value used through- 
out the rest of the present work, Mmm ~ 10~® Mq. 



concentrated in halos with mass between 1.8 x lQ~^Mmr 
(corresponding to the mass resolution in the Aquarius 
simulation) and 10~^Mi,i,.. 

It is convenient to recast the above distribution in 
the form ^(M,r) = N,i'^{M)^{r), where 



M^i„ dM Jo 'tin 



(r), 
This implies 

the definition of a subhalo mass range: while M^ax is 
usually fixed at - IQ-'^Mmw ^ 10^" Mq, Mmin de- 
pends on the nature of dark matter. For supersym- 
metric neutralinos, for example, it varies in the inter- 
val 10"i2 - 10"^ Mq [20., but it may also be much 
larger or smaller. We choose to fix Mmin = 10~^ M© 
which is a typical value for WIMPs [H]. As will be 
shown in sections IV and [V] the positron and antiproton 
clumpy fluxes scale with the quantity Nci{S,)m, where 
(Oai = '^{M)^{M,rQ). Thus, for a given 

subhalo population (either Via Lactea II or Aquarius) a 
different value of Af,„i„ may be accounted for by scaling 
the e'^ and p clumpy contributions according to figure ^ 
Other relevant quantities are the total mass in clumps 
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the local clump fraction 



^^^(r©) and the total clump fraction /*°* = 
Table I displays these and other parameters for 
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bot& Via Lactea II and Aquarius simulations. 



III. PARTICLE PHYSICS INPUT 

The anomalous positron fraction reported by 
PAMELA and the electron plus positron excess 





Via Lactea 11 Aquarius 


Rvir [kpcj 
M^ir [Mq] 

Ts [kpc] 

Ps [10'' Mq kpC-=*] 

Pq [GeV cm ] 


402 433 
1.93-10" 2.50-10" 
21 20 
3.7 2.4 
0.19 0.48 


Or 

Cl 
02 

aci 

Cec2 


0.286 0.237 
119.75 232.15 
-85.16 -181.74 
0.012 0.0146 
0.026 0.008 


[M-' kpc-^] 
Ra [kpc] 

(re) [kpc--'] 
{Om [kpc«] 

M*r*(< R^ir) [Mq] 
fcl (< Rvir) 

fe 


1.7 x 10" 25.86 
21 199 
2.79-10" 1.17-10" 
3.20-10"'^ 8.47-10"** 
3.45-10"" 7.19-10"" 
1.05-10" 4.33-10" 

0.54 0.18 
6.3-10"^ 2.7-10"^ 



TABLE I: Parameters fixing the characteristics of the dark 
matter distribution as deduced from Via Lactea II and Aquar- 
ius results. The three blocks refer to parameters of 1) the 
smooth galactic halo, 2) the concentration parameter and 3) 
the galactic subhalo population. Notice that we are setting 



= 10"^ Mq to compute N^i, (Om, Mlf , /q and 



claimed by ATIC have prompted the interest of the 
particle physics community, and motivated the quest 
for DM models leading to enhanced DM fiuxes. Among 
them, strong emphasis has been put on the Sommerfeld 
enhancement |22j . arising from the presence of an 
attractive potential, that for low relative velocities leads 
to a peculiar behaviour of the annihilation cross section 

O'annV (X 1/v down tO a given Vsat, below which OannV 

saturates. In this scheme, one can have dark matter par- 
ticles that at chemical decoupling — when v ^ 0(0. Ic) 
— presented the appropriate annihilation cross section of 
weak strength, {(Jannv)o 0(10"'^^) cni^s"^, and today, 
in the galactic halo, have a much higher GannV since the 
local velocity dispersion is Z?© = Uq/c ~ 5 • lO""*. The 
subscript in UannV denotes the value of the annihila- 
tion cross section without Sommerfeld corrections. The 
Sommerfeld effect typically leads to small corrections of 
the annihilation cross section at decoupling [23 , while 
boosting significantly local anti-matter fiuxes. Notice 
that the relic abundance of a Sommerfeld-enhanced DM 
particle may be reduced w.r.t. the standard case by a 
factor of order unity: indeed, in Ref. [23 the authors 
have computed the thermal relic abundance of wino-like 
neutralinos including Sommerfeld corrections and found 
a reduction of about 50% compared to the standard 
calculation. 

The Sommerfeld enhancement is rather model- 
dependent and thus we choose a few examples in the 
literature, shown in table [IT] These do not cover all the 
possibilities but are meant to be representative bench- 
marks. We consider specific implementations of Arkani- 
Hamed et al. [14 and Nomura & Thaler [15] models 
as leptophilic-like candidates. In the former, the dark 
matter particle annihilates into pairs of scalars or vec- 
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label 


Ref. 


moM/TeV 


m</,,s/GeV 


(cra„„«)o/(10-^''cm^s-^) 


annihilation channel 


S{vq) 


Smax 


□ AH700 




14 




0.70 


10 (fh) 


3 


fhfh' fh — ^ (^~^ p 


43 


762 


• NTl 




15 




1.00 


34.0 (s) 


3 


SCI] s — > 97%fici, S%bb] a — > 


100 


100 


o NT2 




15 




1.20 


5.60 (s) 


3 


sa; s 95%aa, 5%rr ; a ^ A*"*"/^" 


100 


100 


* fifi 




26 




1.60 




3 




1100 


1100 


M TT 




27 




2.00 




3 


r+r- 


1000 


1000 


X MDM3 




16 




2.70 




~ 1 




273 


273 


® MDM5 




16 




9.60 




~ 1 


WW,ZZ 


1210 


1210 



TABLE II: Properties of DM candidates recently proposed in the literature and presenting Sommerfeld-enhanced cross sections. 
The label of each model is preceded by the corresponding symbol to be used in the plots of the next sections. 



tor bosons cj) which then decay into niuons or electrons. 
We use = 100 MeV (that means 4> decays entirely 
into e+e") and put a = A^/(47r) = 0.01. As for No- 
mura & Thaler models, dark matter annihilates into a 
scalar s and an axion a. We implement the two bench- 
marks of |l5j where s decays mainly into a pair of ax- 
ions but has small branching ratios to bb and tt — see 
table [nj The axion a is assumed to decay entirely into 
muons. Placing Nomura & Thaler models on figure 6 
of [13], one sees that the Sommerfeld enhancement is al- 
ready saturated at w = w© and reads ~ 100 — 300. We 
set S{vq) = Smax = 100 and note that taking a differ- 
ent value is equivalent to rescale {(Jannv)^ since we are 
lying in the saturation regime. Lastly, inspired by the 
recently published electron plus positron spectrum from 
Fermi p3] and HESS [25j, the authors of Ref. [26] pro- 
pose an 1.6 TeV particle annihilating into fj,'^fj,~ , which 
fits well Fermi, HESS and PAMELA data given an en- 
hancement of 1100. Even though such enhancement is 
not necessarily due to Sommerfeld corrections, we con- 
sider this candidate setting S{vq) — S^ax — 1100. Sim- 
ilarly, we also analyse a 2 TeV particle annihilating into 

T+r- with S{vq) = Smax = 1000 pjL 

For hadrophilic candidates we adopt the case of min- 
imal dark matter jl6j, namely the fermion triplet and 
quintuplet. All examples displayed in table [H] are Majo- 
rana fermions. 



Essential ingredients to proceed further are the energy 
spectra dN/dE of positrons, antiprotons and photons 
produced in dark matter annihilations. Arkani-Hamed et 
al. models feature particles that annihilate in an 1-step 
cascade into e~^e~; the relevant formulae for the positron 
spectrum are given in Appendix A of [28i . The Nomura 
& Thaler cases considered here annihilate mainly in an 
1.5-step cascade into ^ which is basically half an 1- 

step cascade and half a 2-step cascade — and the referred 
appendix gives the necessary expressions for dN^+ / dE^+ . 
The corrections due to the branching ratios into bb or 
TT are introduced in the following way: from DarkSUSY 
code version 5.0.4 [29j we take the positron spectra for the 
mentioned channels and ttidm = 10 GeV (corresponding 
to nis — 20 GeV). Then, following [28], we convolute 
them to give the dNf^+ /dE^+ in the centre of mass of 
the DM annihilations. The same procedure is applied for 
dNp/dTp in the bb case. With the small branching ratios 



into TT and bb presented in table |nj the positron spectra 
obtained in this way are very similar to the ones obtained 
in a pure 1.5-step cascade. However, an important dif- 
ference is a non-zero yield of antiprotons (that will turn 
out to be small) and possibly significant 7-ray produc- 
tion. Lastly, energy spectra from minimal dark matter 
annihilations are given in reference [16| for x > 10^*. 
We implement the e+ and p spectrum of [T6 down to 
X — 10^^ and a fiat dN/dE is assumed below that. This 
feature does not affect our conclusions. Indeed, the main 
results will depend on dN^+ /dE^+ for energies above 50 
GeV which are within the range of validity of all spectra 
used. Radio fluxes will depend on the number of elec- 
trons and positrons above ^ 0.1 — 10 GeV; for minimal 
dark matter — where we have implemented a flat dN/dE 
below a given x — the number of will then be under- 
estimated and our radio fluxes for those candidates are 
lower bounds. 



IV. POSITRONS 

Unlike neutral particles, positrons produced in the 
Milky Way undergo different processes that change their 
direction and energy while crossing the galactic medium. 
The galactic magnetic fields, for instance, are responsi- 
ble for defiection and, due to their (poorly known) inho- 
mogeneities, the evolution of a positron can be treated 
as a random walk with a certain diffusion coefficient 
. Other important phenomena are energy losses 
through inverse Compton scattering off the cosmic mi- 
crowave background and starlight and synchrotron emis- 
sion, which proceed at a space-independent rate 6(i?g+) ~ 
El+/{GeV ■ te) with te ^ 10^^ s [3S I3I]- Neglecting 
galactic convective winds and diffusive reacceleration, the 

number density per unit energy ne+(t,x, i?g+) = ^y^"*^ 

follows the diffusion equation [30j [31] 



dt 



-K,+ {E,+ )V^n,+ 



d 



(6(^;,+ )n,+ ) = Q,+ (x, £;,+ ), (2) 



dE,^ 



and we are interested in positrons from annihilations of 
dark matter particles with mass uiem and density pdm 
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L [kpc] Kq 


[kpcVMyr] 


5 


Vc [km/s] 


M2 


1 


0.00595 


0.55 




MIN 


1 


0.0016 


0.85 


13.5 


MED 


4 


0.0112 


0.70 


12.0 


MAX 


15 


0.0765 


0.46 


5.0 



TABLE III: Sets of propagation param eters yi elding maximal, 
mean and minimal anti-matter fluxes I311I32I. 



corresponding to the source term 

' PDM{y^y 



Qe+(x, £'e+) 



moM 



k 



(3) 

where the 1/2 factor is vahd for Majorana self- 
annihilating fermions. The fluxes induced by Dirac 
fermions, bosons or other particles may be obtained from 
our results by a simple rescaling. 

The standard approach to solve (pi for is to as- 
sume steady state conditions (i.e. dn^+/dt = 0) and 
adopt a cylindrical diffusion halo with radius Rgai = 20 
kpc and a half-thickness L inside which the diffusion co- 
efficient is supposed to be space-independent [30l l3T] . 
K^+{E^+) ~ Ko{E^+/GeYf. The half-thickness L ex- 
tends much further than the half-thickness of the galac- 
tic disk ft, ~ 0.1 kpc and 71^+ vanishes at the cylinder 
boundaries since the particles escape to the intergalac- 
tic medium. The propagation model is defined by the 
set of parameters (L, Kq, 6) which turn out to be loosely 
constrained by cosmic ray data, namely B/C measure- 
ments. Following [511 132] we use the sets of parameters 
labelled M2, MIN, MED and MAX in table [HT] that are 
likely to reflect the propagation uncertainty on dark mat- 
ter induced anti-matter fluxes. M2 (MIN) is the set that 
minimises the positron (antiproton) flux. The value of 
the galactic wind speed will be used in the antiproton 
analysis while being neglected here. 

Once the steady state solution is found, the flux of 
positrons is given by 0e+(x, £'e+) = -j^nf.+ (x, Ef.+ ). We 
disregard solar modulation, since it is unimportant for 
multi-GeV positrons. Following [19j [30], the positron 
flux at Earth due to the smooth dark matter component 
of the Milky Way is 
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where the sun is at {xq^uq, zq) — (8,0,0) kpc, Vf.+ /c 



[l-ml/E 



'2^^/'and/,<,(^;s) = E. K 



. dN" 



-{Es). 



\d{E^Es) is the positron diffusion length from a source 
energy Es down to a detection energy E < Es and reads 



Xd{E, Es) 



AKqTe 




5-1 



fEs_\ 
\GeYj 



by 



/|^(A£)) is the dimensionless halo function and is given 



4™(Ad) = / d'x 
Jdz 



Psm (^) 



G^(x,Az,), (5) 



where DZ stands for the cylindrical diffusive zone and 
Gq^ is the Green function evaluated at the solar neigh- 
bourhood: 



Gq^(x, Xd) 



{x - Xq)^ + (y - Vq^ 



with Gfjj given in [19^ (and references therein) for the 
limiting cases L > Xu and L < X^. 

The contribution from one single clump is very similar 
to the smooth one replacing psm with pd in equation (|5| . 
However, we are interested in the signal from a popula- 
tion of subhalos distributed throughout the galaxy in a 



certain range of masses, say 



Ncl 



dPu dPv 



Con- 



dMdV — ^"cl 

sidering every clump a point source and given the local 
character of the Green function, the mean positron flux 
from the clumpy dark matter component in the galaxy is 

m 



i^^JiE) 



Ve+ 1 1 / PG 



An b{E) 2 \mDM 



Ncl (OmX 



dEsflJEs) {G% )v {Xd{E,Es)) , (6) 



where 



{G'q)v{Xd)^ [ d3xG^^(x,Az,) 
Jdz 



dPy 
dV 



(x) 



and (Om was introduced in section |n] Equation (|6| is 
valid if the density profile of the clump does not depend 
on its position within the Milky Way. Therefore, we set 
C20o{M) = C2oo{M,rQ) which is anyway reasonable for 
our analysis since multi-GeV positrons detected at the 
Earth travelled at most a few kpc [30] . 

Next we include Sommerfeld corrections. The smooth 
contribution (f/^^ ^^^^ will be boosted by S{vq) given the 
local origin of high-energy positrons. As far as clumps 
are concerned, we assume that the whole population of 
subhalos presents velocity dispersions below v^at which 
means the clumpy contribution (0^+ will be roughly 
rescaled by Smax- Such simplification is conservative in 
the sense that we maximise the contribution of substruc- 
tures — indeed, clumps with masses close to Mmax = 
10"'^° M0 may not be in the saturation regime, but that 
would lead to an enhancement smaller than Smax- In this 
framework and following [19J, the total positron flux at 
Earth for a specific dark matter candidate and a certain 
propagation model is 



cj,,^{E) = {l-fQYS{vQ)4>\ 



+ ,smiE) 



Srt 



){E). 
(7) 
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FIG. 2: (Color online) The Sommerfeld-enhanced fluxes (1 — 

frs))'^ S{vQ)cj>'^rn and Smax{4>ci) of positrons and antiprotons for 
the MDM3 candidate. Solid red (dashed blue) lines refer to 
Via Lactea II (Aquarius) density profiles. The thick (thin) 
curves represent the smooth (clumpy) contribution. The 
MED propagation model is adopted and Mmin ~ 10~® M©. 



FIG. 3: (Color online) The smooth quantity (1 — 
fQYS{vQ)4)1^ for the MDM3 candidate using the M2 (dot- 
dashed), MIN (dotted), MED (solid) and MAX (dashed) 
propagation models. For clarity, the clumpy component is 
omitted and just the results for Via Lactea II are plotted. 



As an example we show in figure [2] the quantities 
(1 - fQYS{vQ)(lP+^^^ and SmaA€-^J for the MDM3 
candidate presented in table |TTj Both Via Lactea II and 
Aquarius parameters are used and the MED propaga- 
tion set is assumed. On the one hand, the smooth con- 
tribution with Aquarius is larger because its local DM 
density is higher than in Via Lactea II — check table 
|T] On the other hand, the Via Lactea II simulation pre- 
dicts (through extrapolation oi cP N sh / dM dV) more low- 
mass clumps than Aquarius and thus the corresponding 
clumpy contribution is more significant. In fact, the ratio 
between the two clumpy fluxes is simple to understand if 



one computes the ratio 



Aquarius 



0.14. 



Via Lactea II 



Note that, as stated in section |TT] we are fixing Mj^in — 
10~^ M©; clumpy fluxes scale according to figure [l] 

Figure [3] shows the effect of varying the propagation 
parameters for the candidate considered above, for the 
Via Lactea II case without clumps. 
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V. ANTIPROTONS 

As in the previous section, we follow [19] to model 
the flux of antiprotons at the Earth. Differently from 
positrons, antiprotons do not lose much energy by syn- 
chrotron or inverse Compton since rup ^ mg. In- 
stead, they may be swept away by galactic winds, as- 
sumed to be constant and perpendicular to the disk: 
Vc(yi) = sgn{z)Vcez- Furthermore, annihilations pp are 
responsible for the disappearance of primary antipro- 
tons. These annihilations take place essentially along 
the galactic plane where the interstellar medium is con- 
centrated, and therefore the diffusion equation for an- 
tiprotons contains a term —2h5£){z)T'^^-jip with r^fj„ — 

riHe 

and 



661(H-0.0115rp 



36T, 



-0.5 



-0.774 
P 



0.948T? oi5i) Tp < 15.5 
Tp > 15.5 ' 



where ct^J^„ is in mbarn, Tp ~ Ep — nip is the antiproton 
kinetic energy and is in units of GeV in the above for- 



mula. The diffusion coefficient is Kp{Tp) = K^Pp {^^) 

with (3p ={l- (7^^) and Pp = (Tl + 2mprp)i/2. 

All in all, the steady-state solution of the antiproton dif- 
fusion equation obeys 

d 

-Kp{Tp)^^np+— {sgn{z)V^np) = 

Qp{^, Tp) - 2hSD{z)TZ^{Tp)np , 

with Qp of form analogous to ([3| . 

Once again, 0p(x,Tp) = ^np(x,Tp). Because we are 
interested in high-energy antiprotons we have neglected 
solar modulation and reacceleration effects. Similarly to 
the positron case, one has 



with 
and 

«c/)(r) 

where 



Vn 1 



477 2 KniDM 



DZ 



Pe 



fUT)lUT), (8) 



GUx,T,L,Ko,6,Ve), 



vp_l f P& 

477 2 [niDM 



Nci iOMfUT) {Gl)v{T), 
(9) 



(G|)y(r)= / d^KG%{x,T,L,Ko,5,V,)^{x) . 
Jdz 

The Green function for antiprotons Gq is given in |19j . 
At last, the total antiproton flux is 

(bpiT) = (1 - S{vq) 4>lsn.{T) + Sraa. i^JiT) . 

(10) 
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FIG. 4: (Color online) Total positron and antiproton fluxes 
at 50 GeV for the models in table Ull and for both Via Lactea 
II and Aquarius parameters. The MED propagation is used 
and Mmin ~ 10~^ M©. The solid lines indicate the fluxes 
deduced from PAMELA data ass umi ng ele ctron and proton 
fluxes at 50 GeV — see equations (111, ( 12 1 — and the corre- 
sponding excluded regions are shown as snaded. Candidates 
lying above the diagonal line cannot be rescaled to explain 
PAMELA positron excess without overproducing antiprotons. 



Figures |2] and [3] show the antiproton fluxes for the min- 
imal dark matter triplet considered at the end of the pre- 
vious section. 

For clarity in visualising the e+ and p yields from the 
candidates in table |TT] we choose to plot fluxes at 50 GeV 
which is well inside the energy range where PAMELA 
detected the positron excess and collected antiprotons. 
Moreover, both positrons and antiprotons of such en- 
ergy are not signiflcantly affected by solar modulation 
or reacceleration effects. At 50 GeV the PAMELA data 
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indicates [TJ [7] (^g+/(0g- + ~ 0.1 and 4'p/4'p ^ 

0.7 — 2.2 • 10~^ that we translate using figures 9 and 10 
in [3U] and equation 1 in [33] into the fluxes: 



(50 GeV) 



ip(50 GeV) = 



1.0 • 10-^ for hard e" 
4.7 • 10-9 for soft e- 



and (11) 



7.9 
2.5 



10 
10 



"9 for (/)p/0p 
9 for 0p/?;ip 



2.2 ■ 
0.7' 



10- 
10- 



(12) 

in units of GeV ^cm •^s ^sr ^. 

Our aim in the present work is not to perform a fit- 
ting procedure or a likelihood analysis to PAMELA (or 
ATIC) data, but rather to investigate which DM parti- 
cles are able to produce positron fiuxes near the above- 
stated values. Figure |4] shows the total positron and 
antiproton fiuxes at 50 GeV for the models in table |ll] 
and assuming the MED propagation configuration and 
^min = 10"-® M0. We see immediately that some mod- 
els violate the positron and/or antiproton data. How- 
ever, all candidates may present {<Jannv)o a few times 

TT] which accounts 
ITsince a change of 



more or less than presented in table 
for a scaling along diagonals in figure 
(CTanni')o modifics equally the positron and the antipro- 
ton fiuxes. Conservatively, we are interested in knowing 
which points in figure |4] may be rescaled to touch the left 
vertical line without being above the upper horizontal 
line. In other words, we wish to pin down the particles 
that can meet the PAMELA positron excess without vi- 
olating antiproton bounds. Candidates lying above the 
diagonal line in figure |4] cannot. For the DM distribution 
suggested by Via Lactea II the model labelled MDM3 is 
disfavoured. The situation for Aquarius is similar, but 
more constraining. Lastly, leptophilic candidates pro- 
duce no antiprotons and automatically pass the p test; 
they are plotted in figure [4] for completeness and with an 
artificial (t>p. The exception is NTl model that features 
a non-zero BR(s — > hh), even though the corresponding 
antiproton fiux is rather low. 

As mentioned above, figure |4|refers to the MED propa- 
gation parameters. From figure|3|one sees that at 50 GeV 
the positron fiux is not very sensitive to the propagation 
parameters, while the antiproton fiux at 50 GeV may be 
roughly one order of magnitude above or below the fiux 
computed with MED. Since the antiproton bound will 
turn out to be the less constraining one, we can safely 
stick to the mean propagation in presenting our main 
results. 

Next we proceed with constraints coming from 7-ray 
observations. 



VI. GAMMA-RAYS 

The expected 7-ray fiux from DM annihilation is pro- 
portional to the line of sight integral (LOS) of the DM 
density square — which tells us how many annihi- 
lations we have in the cone of view defined by the ex- 



perimental angular resolution — and to the annihilation 
cross-section — which gives us the yields of photons ob- 
tained in one annihilation: 



1 <T„ 



/ 



pUm,r) 



dV . 



(13) 

For each halo along our line of sight, M is the mass 
and d the distance from the observer. dN^/dE^ is the 
differential spectrum per annihilation of prompt photons 
relative to the final state /, with branching ratio 
Bf. The volume integral refers to the line of sight and 
is defined by the angular resolution of the instrument 9 
and by the direction ij: of observation. Inside each halo, 
R is the distance from the centre. py^{M, R, c{M, R)) 
is the DM density profile with concentration parameter 
c{M,r). 

When including the Sommerfeld enhancement of the 
local annihilation cross section, the term UannV in eq. 



( 13 1 is replaced by the velocity-dependent expression 
{crannv)oS{R, M) . The enhancement S depends on the 
halo mass fixing the average velocity dispersion, and from 
the radial coordinate inside the halo, which takes into ac- 
count the features of the velocity dispersion curve that 
has lower values closer to the centre of the galaxy hosted 
by the DM halo - see e.g. [57] . 

The enhancement will be included in the line of sight 
integral, which transforms into 



LOS 



S{M,R)pI{M, c, R) 

d2 



dV . 



(14) 



In the next paragraphs we will compute the prompt 7- 
rays coming from DM annihilation in 1) the smooth halo 
of our galaxy, 2) the substructures of our galaxy and 3) 
the extragalactic halos and subhalos. 



A. Annihilations in the smooth MW halo 

The astrophysical contribution to the prompt 7-ray 
emission from the smooth component of a DM halo can 
be rewritten as the volume integral 



$s(M,d,'(/', Ar2) oc 



d<pded\ 



S{M, R)pI{M, R) 
A2 



where A is the line-of-sight coordinate, AO, the solid angle 
corresponding to the angular resolution 9 of the instru- 
ment, and "0 the angle of view from the GC. Both in 
the case of ACTs and of Fermi (E > 1 GeV) we will use 
9 = 0.1° and Af] - 10"^ sr. 

The HESS telescope has observed the GC source in 
2003 and 2004, measuring an integrated fiux above 160 
GeVof $(> leOGeV) = 1.89 x 10-"7cm-2s-i [38]. In 
figure [5] we show the result of our computation for the 
particle physics models of table |TT| in the cases where 
the Milky -Way halo is described by either the VL2 or 
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the Aquarius models. Also shown is a table with the 
expected flux above 160 GeV at the GC, for direct com- 
parison with the HESS limit. Considering {aannv)o — 



B. Annihilations in the subhalos of the MW 



We now populate the halo of our galaxy with subhalos 
with masses as small as 10~^ Mq. 

The LOS contribution of such a population of substruc- 
tures to the annihilation signal can be written as 



<l>^s"^(A4,d,V,Af])cx 



dM. 



sub 



M, 



dc 



AO 



dX[p,h{Mh,Msub,R)P{c)^T"'] 



(15) 



where the contribution from 



$^°'°(M, 5(M, R),c{M, R),d, An) 



each subhalo 

yivi ^ ijyivi ^ ih) ^ t_^yivi ^ iLj ^ u,^ Ljj ^ L-^iL j is convolved 
with its distribution function inside the galaxy psh- Here 
R is the radial coordinate with respect to the centre of 
the host galaxy. -P(c) is the lognormal distribution of 
the concentration parameter with dispersion = 0.14 
[39l and mean value c: 



P(c,c) 



1 



Mn(c)-ln(c) ' 



(16) 



The satellites inside a galaxy suffer from external tidal 
stripping due to the interaction with the gravitational 
potential of the galaxy itself. To account for gravitational 
tides, we follow [ID] and assume that all the mass beyond 
the subhalo tidal radius is lost in a single orbit without 
affecting its central density profile. The tidal radius is 
defined as the distance from the subhalo centre at which 
the tidal forces of the host potential equal the self-gravity 
of the subhalo. In the Roche limit, it is expressed as: 



\2Mh{< r) 



1/3 



(17) 



where r is the distance from the halo centre, Msub the 
subhalo mass and Mh{< r) the host halo mass enclosed 
in a sphere of radius r. The integral along the line of 
sight will be different from zero only in the interval [d — 
rtid,d + rtid]- 

The estimate of the 7-ray flux from substructures is 
obtained through the numerical integration of Eq. ( |15| . 
The result of our computations is shown in flgure[6] where 
the expected flux from the subhalos is compared with the 
EGRET diffuse emission (galactic -I- extragalactic) back- 
ground. Although the EGRET flux has been rescaled 
by ~15% in the light of the new measurement from the 
Fermi telescope which does not confirm the GeV bump, 
the diffuse emission measurements turn out to be less 
constraining than the GC one. 



VL2 



HESS GC 




5.0 10.0 

if/ [deg] 




5.0 10.0 
(A [deg] 



<f>Mw{lp = 0°) 



label 


Via Lactea 11 


Aquarius 


AH700 

NTl 

NT2 

nn 

TT 


1.15 ■ 10"'^ 
1.33 ■ 10"" 
2.39-10"" 
8.67-10"^^ 
4.47- 10"" 


3.38- 10"'^ 
4.16 ■ 10"" 
7.47- 10"" 
2.70 ■ 10"" 
1.39 ■ 10""' 


MDM3 
MDM5 


2.01 ■ 10"'^ 
2.47- 10"^^ 


6.25 ■ 10"'^ 
7.69- 10"^^ 



FIG. 5: (Color online) The 7-ray flux above 160 GeV as a 
function of the angle with respect to the GC. The legend in 
the figures is ordered according to the values of the curves at 
ijj = 0.1°. The HESS measurement towards the GC is shown 
by the horizontal thick line and the table shows the fluxes for 
iP = 0°. 



10 



^ 01 



^ 0.01 
•§0.001 



10--* 



10- 



I 10- 



o 
Si 



10" 



^ 0.1 



VL2 




MDM5 








"TTrrr^^.v;-.^^^^ _ TT 


w _ - " 












NT2 







^EGRETiE = lOGeV) = 6 X lO'^ cm-^ g-i. 

The prompt photon flux from DM annihilation in the 
extragalactic host halos can be written as: 
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FIG. 6: (Color online) Number of photons above 3 GeV 
from unresolved subhalos at a galactic longitude 1 = 0° 
for Via Lactea II and Aquarius parameters. Again, we fix 
Mmin = 10"'' M0. Shown in thick solid are the EGRET map 
for diffuse background as well as the smooth approximation 
of reference 141], both scaled to Fermi results. 



C. Annihilations in the extragalactic halos and 
subhalos 



Any direction of observation will also receive a con- 
tribution from all the halos and subhalos at all red- 
shift along the line of sight. We adapt the formalism 
of Ref. [42] to estimate their contribution to the prompt 
7-ray flux in case of Sommerfeld enhancement. We com- 
pute the extragalactic contribution at E = 10 GeV, and 
compare it with the EGRET extragalactic background 
at the same energy, rescaled for a 40% due to the contri- 
bution of those unresolved blazars which have been now 
observed and included in the Fermi catalogue, namely 



dEodn^ dEodtodAdn 



(18) 



where Eq and dto are the energy and the time interval 
over which the photons are detected on earth. dN^ is the 
number of 7-ray photons produced in the infinitesimal 
volume dV at redshift z {dV cx dil{l + z)^'^), in a time 
interval dt [dt = dtQ{l + z)~^) with an energy between E 
and E + dE (E — Eq{1 + z)) and collected by a detector 
with effective area dA. 

dNj is obtained integrating the single halo contribution 
to the photon flux dAf^/dE over the halo mass function 
JMiM,z): 



dN 



-y cx e 



-t{z,Eo) 



{1 + zr 



dM^{M,z) (19) 



dM.y 

~dE 



{E,M,z) 



Following the Press-Schechter formalism ^43^, the halo 
mass function is given by 

dn , , Pcrf^o.m dlogiy 
— (M,z) = ^^./(.)-^^, (20) 

where p^, is the critical density, fio.m is the mass 
density parameter, v = ^^m) ' cr{AI) is the rms density 
fluctuation on the mass scale M and 6sc represents the 
critical density for spherical collapse |1H HH |3S] . 

dAf-y I dE is the number of photons with energy between 
E and E + dE produced in a halo of mass M at redshift 
z, and can be written as 



dAf-y 
~dE 



{E,M,z) = 



dN^{E) M Ay„Pcr^^{z\ 



2 E 

c^{M,z) 



'{21) 



-/2(2:.„.„,c(Af, z)). 



h{x^,^,c{M,z)Y 
In the previous expression, the virial overdensity is |42J: 
IStt^ + 82(r2„(z) - 1) - 39(r2„(z) - 1)2 

(22) 



and the integrals /i and I2 enter the LOS integral 

In{x^,n,x^..) = I 9"x^dx, (23) 



where g{x) = x ^(1 + a;) ^ in the case of an NFW profile 
and g{x) = e^^C^; -i)/" j-^ ^j-^g Einasto case. 

Finally, the absorption coefficient g-'^i.^^Eo)^ 
t{z,Eo) = z/[3.3(i;o/10 GeV)-°-8] accounts for pair 
production due to the interaction of the 7-ray photons 
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with the extra-galactic background light in the optical 
and infrared bands |46j . 

With a little algebra, we get to the final expression for 
the extragalactic DM 7-ray flux: 



C 



2 



(24) 



hyz) dE 



scattering contribution is particularly relevant for lep- 
tophilic models and has been computed in |47l [48] where 
it has been shown that COMPTEL and EGRET extra- 
galactic observations place interesting limits on DM an- 
nihilation cross-sections. These constraints are competi- 
tive with the ones derived here with prompt 7-rays from 
the Galactic Centre. Another effect of the presence of 
inverse Compton photons is the ionisation of the bary- 
onic gas after recombination and thus the decrease of the 
CMB optical depth - see Refs. [49ll50]. 



with 
and 



dM 



v{z,M)f{y{z,M)) 
Ma{M) 



da 



dM 



MAli{z,M) 



(25) 



Aii{z,M)^ dc'P{c{M,z),c') 



, A„i,. /2(2:mi„,c') , ,,3 



3 /?(x_,c') 



(26) 

A'^(z) is a quantity which describes the boost to the 
isotropic 7-ray background due to the existence of 
virialized DM halos. 



To account for the presence of substructures inside the 
extragalactic host halos, we have to apply the substitu- 
tion 

MAl,{M) -> (l-/(M))2MA2,(M) + Ai,,,„,(M) (27) 

where f{M) is the fraction of mass in virialized substruc- 
tures within a host halo of mass M , and can be expressed 
as 



VL2 



1.45 X 10"^[Zo5(M) -f 9.21] , 



f{M)Aq = 0.285Af-"-^[0.63M''-^ - 0.25] . 

More precisely, taking into account the effect of the ra- 
dial dependence of the concentration parameter of sub- 
halos, we will have 



^\l,subi.Z, M) 



dMsub / dc'P{c{M,z),c') (28) 



Airr'^drpshiM, r) 



W /2(a:.„in, c'(Af,r)) 
3 /i2(x„,„,c'(M,r)) 



(c'(M, r))3dc'. 



We obtain our results by solving numerically Eq. (24) 



both for host halos and for subhalos, and combining them 
through Eq. ( 27 1 . The extragalactic bounds turn out to 



be the less constraining. Note that our bounds refer to 
prompt 7-rays only. However, DM annihilations in extra- 
galactic halos and subhalos can also give rise to a 7-ray 
flux by producing high-energy electrons and positrons 
that up-scatter CMB photons. This inverse Compton 



VII. SYNCHROTRON RADIATION 

Synchrotron emission arises from relativistic electrons 
and positrons propagating in the galactic magnetic fleld. 
Since all annihilation channels usually considered pro- 
duce high-energy electrons and positrons, a DM- induced 
synchrotron signal is expected from regions of the galaxy 
where a substantial magnetic field is active and the dark 
matter density is significant — see e.g. [3SJ [SD El] . Let 
us focus on a region towards the galactic centre, small 
enough so that diffusion does not play an important role 
and where the galactic magnetic field is strong enough to 
neglect electron (and positron) energy losses other than 
synchrotron emission. Assuming further that advection 
is negligible as in |36j |52] , equation ^ is easily solved in 
steady-state conditions for positrons and electrons: 



ne±(x,i;) 



v) 2 .N,±{>E) 
—PdmW] 



.(x,£;) ' 



(29) 



where A^e±(> E) = 



^rnoM ^^(^E') is the number 
of electrons plus positrons per DM annihilation above a 
given energy E, bsyn{^,E) ~ e'^B^{ii)E^ /{9T:m^) and B 
is the galactic magnetic field. Now, each electron (or 
positron) emits synchrotron radiation according to the 
power spectrum 15^ 



dWs^ 



dv 



:(x,£;) 



V3e3B(x) 87r 



9^3 



syn 



0.29 



(30) 

being i'syn{x,E) = 3eB{x)E'^ /{Anml) the synchrotron 
frequency. The next step is to convolute equations ( 29 1 
and ( 30 ) in the volume of observation Vobs , in order to 



give the total synchrotron power emitted by the distri- 
bution of DM-induced electrons and positrons: 



dW,^, 



j dV J dEn,±{x,E)v^^^{x,E) 
dV PDMWEp{x,i^) ^,(31) 



diy 



where Ep{x,v) = ^ 3 0.29eB(x) 

It was pointed out in |52j that low- frequency radio con- 
straints do not depend much on the magnetic field pro- 
file adopted. Following that work, we choose to imple- 
ment a constant B = 7.2 mC inside the accretion region 
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r < Race = 0.04 pC, B (X for Race < r < 84.5Raee 

and B = for r> 84. 5i?acc- 

In order to place constraints on "^"2 , we consider 

the three configurations studied in [53]: a cone of half- 
aperture 4" pointed at the GC and v = 0.408 GHz (case 

1) , a region with angles from the GC between 5' and 
10' and I' = 0.327 GHz (case 2), and finally a cone of 
half-aperture 13.5' pointed at the GC and v = 0.327 
GHz (case 3). We use the measured fiuxes quoted in 
[53J. It turns out that, independently of the annihilation 
channel, case 1 gives the most stringent bounds when 
using the NFW smooth profile suggested by Via Lactea II 
simulation. For the Aquarius simulation and its Einasto 
smooth profile, case 2 is the most constraining one. 

In principle, for DM particles with Sommerfeld- 
enhanced cross sections, a full calculation of the syn- 
chrotron emission should include the enhancement S{v) 
inside the integral in expression (31 1. Nevertheless, 
for our present proposes it suffices to put GannV ~ 
S{vQ){<Jannv)Q siucc the signal comes mainly from re- 
gions where u ~ uq — recall that v ~ at r = Race = 
0.04 pc. In case 1, for instance, one is looking into a re- 
gion of size ~ 0.16 pc around the GC. The region defined 
in case 2 encompasses distances of ~ 10 — 20 pc from 
the GC where v < Vq and S{v) > S{vq); hence, in this 
case, considering UannV ~ S{vQ){(Jannv)Q yields actually 
a lower bound on the radio fiux. 

The results for the radio flux (case 1 for Via Lactea 
II and case 2 for Aquarius) are presented in figure [t] 
plotted against the corresponding positron fluxes at 50 
GeV. In the case of Via Lactea II, we can see that most 
candidates seem to be at odds with radio observations 
even when a rescaling of (cra„„w)o is applied to meet the 
positron excess. The situation for Aquarius (using case 

2) is similar. 



VIII. CONCLUSIONS 

Table |IV| summarizes our main results. There, we 
present for each model under consideration the max- 
imum (tTa„„w)o allowed by the antiproton bound, the 
HESS measurement from the GC and radio observations. 
For antiprotons we conservatively use the largest value in 
equation (12 1. The most constraining of the three limits, 
i.e. the one yielding a minimal {a an 



iv)a,-,nax, is displayed 
in bold. We have disregarded here diffuse 7-rays as well 
as radio fluxes in the cases 2 and 3 (1 and 3) for Via 
Lactea II (Aquarius) since they give subdominant con- 
traints. Furthermore, the column shows the value of 
(farm*^)o needed to meet the lowest positron flux in equa- 



tion (111; these numbers arc underlined only if allowed 
by the most stringent bound in bold. Notice that we 
apply this procedure to MDM3 and MDM5 even though 
minimal dark matter is a rather predictive scheme. 

Firstly, we immediately see from table |IV| that radio 
observations are rather constraining w.r.t. antiprotons or 
7-rays for models trying to explain the positron excess. 
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FIG. 7: (Color online) The radio flux for case 1 (2) and Via 
Lactea II (Aquarius) parameters against the positron flux at 
50 GeV already presented in figure|4] The MED propagation 
is used and Mmin = 1O~^M0. The horizontal line indicates 
the measured flux (see |53) for details). Similarly to figure H] 
candidates lying above the diagonal line cannot be rescalea 
to explain PAMELA positron excess without violating radio 
bounds. Shaded regions are excluded by PAMELA (assuming 
the electron fluxes discussed in the text) and radio observa- 
tions. 



In the case of Via Lactea II DM distribution, just one 
candidate survives the studied bounds: AH700. There is 
still some tension with the data though. In fact, in order 
to flt the PAMELA data, this model needs to be rescaled 
down to {(JannV^Vi of ~ 10~^^ cm'^s"^, which means ei- 
ther that such particles would overdose the universe in 
the standard thermal relic scenario, or that they are pre- 
vented to be the dominant dark matter component in the 
universe. As far as Aquarius is concerned, we identify no 
model that can evade all the implemented bounds. 
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FIG. 8: (Color online) Multi-messenger bounds for several DM candidates. The MED propagation is used and Mmin = 
10~^ Mq. The e+, p, 7-ray, radio (cases 1 and 2) axes are normalised to (t>^+{50 GeV) = 7.35 • 10"*' GeV"^ cm ^s ^sr ^, 
0p(5O GeV) = 5.2 • 10"^ GeV-^cm-^s-^srli (t)MwUb = 0°) = 1.89 ■ 10""cm-2s-\ 0.05 Jy and 121 Jy, respectively. The 
boxes encompass the values in equations (111 and (12 I for positrons and antiprotons, and a 20% uncertainty on top of the 
HESS measurement for 7-rays. Notice that changing the value of {oannv)o leads to an overall scaling of the polygon. 
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(cranuW)o,maa:/(10 "'"cm'' 








Via Lactea 11 


Aquarius 


label 


e+ 


P 


7 GC 


radio (1) 


e+ 


P 


7 GC 


radio (2) 


AH700 


U.34 




49 


2.0 


1.2 




17 


0.76 


NTl 


4.4 


840 


427 


1.2 


2.2 


200 


136 


0.51 


NT2 


6.0 




238 


1.8 


2.9 




76 


0.74 


MA* 


1.8 




6.5 


0.80 


0.81 




2.1 


0.31 


TT 


3.1 




1.3 


1.1 


1.5 




0.41 


0.43 


MDM3 


12 


7.9 


9.4 


0.86 


5.7 


1.9 


3.0 


0.54 


MDM5 


9.9 


69 


7.7 


1.7 


4.9 


18 


2.5 


0.82 



TABLE IV: The maximum allowed {(Tannv)o by the antiproton bound, the 7-ray measurement from the GC and radio observa- 
tions. The MED propagation is used and Mmin = 10~^ Mq. The bold values represent the most constraining channel and the 
columns display the {aannv)o needed to meet the positron excess. Underlined values manage to explain the positron excess 
while being allowed by our multi-messenger scheme of constraints. 



In figure [8] we introduce a new method to visualise the 
multi-messenger constraints. We place each of the chan- 
nels in a semi-axis and normalise it to the experimental 
limits. Since we are interested in the viability of the DM 
explanation of the positron excess, we only show models 
that are able to reproduce the observed PAMELA flux, 
and thus cross the "up" axis at 1. Configurations exceed- 
ing the boxes on other axes violate observational bounds 
and are therefore ruled out. Configurations not crossing 
the boxes are in principle viable, but one has to check 
then whether the cross section allows to achieve the ap- 
propriate relic abundance, whether the "boost-factors" 
are reasonable, and whether the model provides a good 
fit to all PAMELA data. We stress that once the model 
is specified, the shape (angles and number of vertices) is 
fixed, and changing the cross-section corresponds to in- 
creasing or decreasing the overall size of the polygon. 
Note for instance that leptophilic configurations (first 
and second rows of figure [8| have a different shape with 
respect to hadrophilic ones (third row). 

To sum up, we have analysed the possibility that the 
models in table |ll] recently suggested in the literature, 
may explain the PAMELA positron excess without 
violating bounds in the antiproton, 7-ray and radio 
channels. It turns out that — even considering both 



substructure and Sommerfeld enhancement — the candi- 
dates that provide a good fit to positron data, inevitably 
overproduce antiprotons, gamma-rays or radio emission. 
Our conclusions hold for the DM distributions from 
Via Lactea II and Aquarius simulations, the MED 
propagation model and M,nin — 10^^ Mq. As discussed 
in section |V] modifying the propagation parameters 
does not change our main results. This is because the 
antiproton bound — the most sensitive to propagation 
— is subdominant. Thus, a different propagation model 
is not sufficient to reconcile the studied candidates with 
the observational constraints. A non-standard DM 
profile, and a non-standard magnetic field profile at 
the Galactic Centre, can in principle make theoretical 
models compatible with observations, at the expenses of 
introducing new ad-hoc hypotheses on these quantities. 
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